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FOREWORD 


The  Professional  Development  Center  is  the  primary  source  of  young 
civilian  engineers  and  scientists  for  the  Naval  Air  Systems  Command. 

At  one  point  in  their  training  program,  they  undertake  an  original  Special 
Project  as  part  of  the  requirements  for  an  accelerated  promotion.  Some 
of  the  reports  on  these  special  projects  have  been  both  interesting  and 
informative,  and  deserve  somewhat  wider  distribution.  The  results  pre¬ 
sented  herein  are  not  intended  to  reflect  official  US  Navy  policy,  nor 
necessarily  even  the  views  of  the  Naval  Air  Systems  Command.  The  results 
of  the  Special  Project  are  presented  herein  because  they  are  interesting, 
and  because  they  may  constitute  a  small  contribution  to  the  literature. 
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ABSTRACT 


Non-combat  aircraft  accident  statistics  indicate  that  a  direct 
relationship  exists  between  the  number  of  accidents  and  accumulated 
flight  hours  or  similarly  between  the  accident  rate  and  accumulated 
flight  hours  for  each  model  of  military  airplane.  This  paper  inves¬ 
tigates  the  feasibility  of  relating  accident  rates  directly  to  the 
total  number  of  past  accidents. 

Based  on  the  pure  birth  process  a  method  for  predicting  aircraft 
accidents  is  presented.  Application  of  this  procedure  to  various  test 
cases  shows  interesting  and  useful  results.  One  definite  conclusion 
that  can  be  drawn  is  that  with  two  or  more  years  of  flight  and  accident 
data,  future  aircraft  accident  rates  can  be  predicted  with  fairly 
reliable  accuracy. 

Future  studies  based  on  these  same  procedures  will  delve  into 
further  relationships  that  may  exist  between  aircraft  characteristics 
and  other  relevant  accident  factors. 
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INTRODUCTION 


United  States  Navy  and  Air  Force  statistics  of  non-combat  air¬ 
craft  accidents  indicate  that  for  each  model  aircraft  some  direct 
relationship  exists  between  total  number  of  accidents  and  accumulated 
flight  hours  or  equivalently  between  accident  rate  and  accumulated 
flight  hours.  Many  studies  have  been  conducted  in  the  past  to  discover 
the  nature  of  this  relationship.  The  present  study  investigates  the 
possibility  of  relating  accident  rates  directly  to  total  number  of  past 
accidents  (instead  of  accumulated  flight  hours).  The  relationship 
between  accident  rate  and  accumulated  flight  hours  thereby  appears  only 
as  an  indirect  consequence  of  the  relationship  which  exists  between 
accident  rate  and  number  of  past  accidents. 

A  specific  method  for  predicting  aircraft  accidents  is  proposed 
based  on  the  pure  birth  process.  A  sample  case  is  set-up  and  run  to 
demonstrate  the  usefulness  of  the  theory  and  the  computer  programs. 
Alternate  approaches  to  the  problem  are  presented  for  comparison  and 
evaluation  purposes.  The  method  used  for  a  specific  case  may  depend 
on  the  trends  demonstrated  in  the  data. 

All  references  in  this  report  to  aircraft  accidents  apply  to  non¬ 


combat  aircraft  accidents  unless  otherwise  stated. 


I.  THE  THEORY  AND  ASSUMPTIONS 


The  U.  S.  Naval  Aviation  Safety  Center's  statistics  for  aircraft 
accidents  are  presented  by  quarters  for  the  years  1954  to  1962  and 
annually  thereafter  up  to  the  present.  For  each  reporting  period,  the 
number  of  flight  hours  or  landings  and  corresponding  number  of  accidents 
under  various  classifications  (by  aircraft  model,  damage  and  injury 
classes,  fleet,  etc.)  are  tabulated  in  the  reports.  The  Safety  Center's 
reports  have  been  undergoing  continual  improvements  and  expanded  cover¬ 
age  over  the  years  so  that  there  are  certain  items  found  in  later  reports 
that  are  missing  in  earlier  reports.  For  purposes  of  this  unclassified 
report,  it  shall  be  assumed  that  the  flight  hours  T'  and  number  N'  of 
accidents  of  a  specific  aircraft  model  for  a  statistically  significant 
number  m'  of  consecutive  report  periods  can  be  extracted  from  the  Safety 
Center's  reports  and  can  be  displayed  as  follows. 


TABLE  1.  Initial  Data 


Report  Period 
Number 

1 

2 


Flight  Hours 


T' 

2 


Number  of 
Accidents 


N' 

2 


m' 


T' 

m 


N' 

m 


The  report  periods  are  not  necessarily  of  the  same  leng.h.  This  section 
discusses  the  basic  statistical  model  assumed  in  this  report  for  the 
analysis  of  such  data. 
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The  cumulative  number  n'  of  accidents  and  flight  hours  t'  are 

defined  for  K  =  0,1, 2,..., m*  by  n'  *  0,  t'  *  0  and  for  0  <  Kim', 

0  0 


ni  -  N!  +  N'  +  ...  +  Ni  , 


ti-T;  +  Ti  +  -  +Ti- 


When  n'  is  plotted  against  t',  the  resulting  points  usually  appear  to 
fall  in  a  neighborhood  of  a  continuous  curve,  (See  Figure  1.)  The 
object  of  this  study  is  to  find  a  method  of  defining  the  underlying 
curve  so  that  conclusions  regarding  accident  rates  may  be  derived  from 
it. 

The  method  of  analysis  employed  in  this  report  requires  that  for 
each  period  j,  j  ■  1,2,3,...,  the  number  of  accidents  be  positive, 

Nj  >  0,  The  data  is  therefore  modified  to  eliminate  any  periods  where 
Nj  ’  *  0  by  the  following  rules.  If  the  number  of  accidents  begins  with 
a  string  of  zeros,  N|  *  *  . . .  “  *  0,  followed  by  N^+1  >  0,  set 


T.  +T.  +  ...  +T.  +  t£+1 


as  the  flight  hours  for  a  combined  first  period  with  accidents. 

If  N^ '  >  0  is  followed  by  a  string  of  zeros,  “  ^j+2  "•••"  **j+K  *  ® 

and  N'  >0,  set 
j+K+1 

T'  +  T'  +  ...  +  T' 
j+1  j+2  j+K+1 


as  the  flight  hours  for  a  combined  period  with  Nj+£+^  accidents.  If 
Nj  >  0  is  followed  by  a  string  of  zeros,  Nj+1  '  8j+2  NJ+K  '  0 
V  is  the  last  entry,  disregard  all  data  after  period  j. 
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FIGURE  1 


INITIAL  ACCUMULATED  DATA  GRAPH 


CUMULATIVE  AIRCRAFT  HOURS 
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Since  the  U.  S.  Navy  statistics  on  accident  rates  are  given  in 
units  of  the  number  of  accidents  per  10,000  flying  hours,  the  data  must 
be  further  modified  so  that  the  flight  hours  are  reduced  to  units  of 
10,000  hours. 


TABLE  2. 

Modified  Data 

Period 

Flight  Hours 
x  10- 4 

Number  of 
Accidents 

Rate 

1 

T1 

N1 

Nj/T 

2 

• 

• 

T2 

• 

• 

N2 

• 

• 

VT 

• 

• 

• 

m 

• 

T 

m 

• 

N 

m 

• 

N/T 
m  ] 

Finally  for  K  »  0,  nQ  *  0  and  tQ  =  0,  and  for  0  <  K  <.  m, 


nK  "  Ni  +  N2  +  •  •  •  +  %, 

tR  -  I,  +I2  +  ...  +Tr 


The  plot  of  tK  versus  n^  (see  Figure  2)  looks  like  the  original  plot  of 

t'  versus  nA  (Fibure  1)  except  that  points  (t!,  n!)  which  show  no  increase 
K  k  j  J 

of  n'  with  time  t!  have  been  eliminated  and  the  accumulated  time  has 
K  J 

been  altered  by  the  factor  10“4. 

Now  that  the  problem  to  be  investigated  is  set  up,  it  is  necessary 
to  explain  the  theory  that  applies  to  this  type  of  situation. 

A  stochastic  process  is  an  indexed  family  of  random  variables  Xt  on 
a  probability  space  with  index  t  ranging  over  a  suitable  parameter  set  T. 
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FIGURE  2 


MODIFIED  ACCUMULATED  DATA  GRAPH 


0  *1*2  fc3  S  C5  C6 

CUMULATIVE  FLIGHT  HOURS  X  10'4 


6 


The  state  space  of  the  process  is  a  set  S  in  which  possible  values  of 
each  Xt  lie.  In  the  particular  case  being  dealt  with  in  this  report, 

S  *  {0,1,2,...}  where  the  integers  0,1,2,...  represent  the  accumulated 
number  of  accidents,  so  the  process  is  called  integer  valued  or  a  dis¬ 
crete  state  process.  If  T  *  [Q,“),  as  in  this  problem,  where  t  is 
interpreted  as  accumulated  time,  then  Xt  is  a  continuous  time  process. 

A  sample  function  of  a  stochastic  process  {Xt,  t  e  T}  is  an  assignment, 
to  every  t  c  T,  of  a  possible  value  of  Xt.  Given  the  value  of  Xt,  such 
that  the  values  of  X8,  s  >  t,  do  not  depend  on  values  of  Xu,  u  <  t,  then 
the  stochastic  process  is  Markovian.  That  is,  a  process  is  Markovian 
if  the  probability  of  any  particular  future  behavior  of  the  process, 
when  its  present  state  Is  known  exactly,  is  not  altered  by  additional 
knowledge  concerning  its  past  behavior: 

Pr  (a  <  Xt  <  b  |  Xt  =  x  ,  Xt  *  x  ,  . . . ,  Xt  *  xn} 

1  1  2  1  n 

■  Pr  {a  <  Xt  <_  b  |  Xtn  ■  xn  } 

where  t  <  t  <,..<  t  <  t.  The  function 
1  2  n 

P(x,s;t,A)  =  Pr  {Xt  eA  j  Xs  *  x},  t>s 

is  called  the  transition  probability  function.  A  Markov  process  has 
stationary  transition  probability  if  P(s,x;  t,A)  is  a  function  only  of 
t  -  s.  For  the  special  case  where  A  is  the  one  point  set  {j}, 

Pij(t)  -  Pr{X(t+u)-j|X(u)-i),  i, j *0,1,2, .. . 

is  the  transition  probability  function  for  t>0  and  is  independent  of 
u  >  0. 
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One  example  of  a  continuous  time,  discrete  state,  Markov  process  is 
the  Poisson  process.  If  the  sample  function  Xfc  counts  the  number  of 
times  a  specified  event  occurs  during  the  time  period  from  zero  to  t, 
then  each  possible  Xt  is  represented  as  a  nondecreasing  step  function. 
The  specified  event  occurs  first  at  time  tJf  then  at  time  t2,  at  time 
t3,  etc.;  so  the  total  number  of  occurrences  of  this  event  increases 
only  in  unit  jumps,  and  Xq  -  0. 

The  postulates  relating  to  the  Poisson  process  are: 

(1)  The  number  of  events  happening  in  two  disjoint  Intervals 
of  time  are  independent.  Suppose  tQ  <  tj  <  t2  <• . .<  tfl,  then  increments 
Xtl  -  Xto>  Xtz  -  Xtl,...,  Xtn  -  are  mutually  independent  random 

variables. 


(2)  Random  variable  XtQft-XtQ  dePends  only  on  t  and  not  on  tg  or 
on  the  value  of  Xt0. 

(3)  Probability  of  at  least  one  event  happening  in  a  time 
period  of  duration  h  is 


p(h)  -  Pr{X(t+h)  -X(t)  -  1  |  X(t)  -  x} 

■  Xh  +  o(h) ,  X  >  0. 


(4)  Probability  of  two  or  more  events  happening  in  time  h  is 
o(h).  This  excludes  the  probability  of  the  simultaneous  occurrence  of 
two  or  more  events.. 

(5)  X(0)  -  0. 

Using  these  five  postulates  it  can  be  proven  that  Xt  has  a  Poisson 
distribution  with  parameter  Xt  for  every  t  as  shown  by  Karlin  in  A  First 
Course  in  Stochastic  Processes  (pp  14-16) . 
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Let  Pj(t)  denote  the  probability  that  exactly  j  events  occur  in 
time  t, 

Pj(t)  -  Pr  (Xt  -  J},  j-0,1,2 . 

Postulate  (4)  can  be  written  in  the  form 

E  P,(h)  -  o(h) 
j-2 

and  clearly 

p(h)  -  Pt(h)  P2(h)  +  ...  . 

Due  to  the  assumption  of  independence  in  Postulate  (1) 


-■£&)-  +  X  . 
h 

Therefore  probability  P0(t)  that  the  event  has  not  happened  during  (0,t) 
satisfies  the  differential  equation 

P0'(t)  -  ~XPq (t) 

whose  solution  is 

PQ(t)  -  ce"At. 

9 


10 


Therefore 


^L~")h~  Pj<t)  +  XPJ.1(t)f 


as  h  ♦  0 


resulting  in 


V  3-1.2.  •• 

which  is  subject  to  the  initial  conditions 


Pj  (0)  -  0, 


j  -  1,2,  .. 


To  solve  this  last  differential  equation  substitute 

Qj(t)  ■  Pj(t)  e**,  j-  0,1,2,. 

into  the  differential  equation  P ^ ’  (t). 

Then 

V(t)  ■  3-1.2... 

where 

Q0(t)  -  1 

and  the  initial  conditions 

Qj(0)  "  °*  5  -  1,2 . 


Solving  Qj'(t)  recursively 

Qx’(t)  ■  A  or  Qj (t)  -  At  +  c 


so  Qi<t)  -  At, 


Therefore 

Pj(t)  -4r  e'Xt  . 

J  j  • 

That  is,  for  every  t,  Xt  follows  a  Poisson  distribution  with  parameter 
Xt. 

The  pure  birth  process  is  a  generalization  of  the  Poisson  process 
in  which  the  chance  of  an  event  occurring  at  a  given  instant  of  time 
depends  on  the  number  of  events  that  have  already  occurred. 

A  Markov  process  satisfying  the  following  set  of  postulates  is 
termed  a  pure  birth  process  as  stated  by  Karlin  in  A  First  Course  in 
Stochastic  Processes: 

(1)  Pr{X(t+h)  -X(t)  -  1  |  X(t)  -  k}  -  Xkh+o1>k(h)  (h-0+) 

(2)  Pr{X(t+h)  -X(t)  -  0  |  X(t)  -  k}  -  l-xkh+o2>k(h) 

(3)  X(0)  -  0 

(4)  Pr{X(t+h)  -X(t)  <  0  |  X(t)  -  k}  -  0,  (k^O). 

Define 

Pj(t)  -  P{X(t)  -  j}. 

A  system  of  differential  equations  satisfied  by  Pn(t)  for  t  >  0 
can  be  derived: 

po'(t)  --\)ro(t>* 

Pj'rt)  -  -XjPjCt)  +  »,_!?.  jCt),  j  i  1 

with  boundary  conditions 
PQ(0)  -  1 
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and 


Pj(0)  -  0, 


j  >  0. 


If  h  >  0,  j  2l  1  and  by  use  of  the  law  of  total  probabilities,  the  Markov 
property,  and  Postulate  (4) 

Pj (t+h)  -  Z  Pk(t)Pr  (X(t4h)  -  j  |  X(t)  -  k} 
k“0 

CO 

-  Pk(t)Pr  {X(t-rti)  -x(t)  -  j-k  I  x(t)  -  k} 

-  I  Pk(t)  Pr  (X(t+h)  -X(t)«  j-k  I  X(t)  -  k}. 
k-0 

For  k  -  0,1,2,. ,  n-2 

Pr  (X(t+h)  -X(t)  -  n-k  |  X(t)  «  k} 

<  Pr  {X(t+h)  -X(t)  ^  2  |  X(t)  -  k} 

■°l,k<h)  +  02,k(h> 

or 

Pr  {X(t+h)  -X(t)  -  j-k  |  X(t)  -  k}  -  o3>j(k(h). 

So 

Pj  (t+h)  -  Pj(t)  [1-Ajh  +  o2  j(h)]  +  Pj^(t)  [Aj^h 

+  +  li  Pk<‘>  °3.J,k<h) 

or 

Pj  (t+h)  -Pj  (t)  -  Pj(t)[-Ajh+o2  J(h)]  +  *j-1(t)lAj_1h 

+  oitj_i(h)]  +  Oj (h) 
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where 


11a  — ii^..  -  0 
h-+0  h 


n-2 

uniformly  in  t  >  0  since  on(h)  Is  bounded  by  the  finite  sum  £  o,  ,  k(h) 

k»0 

which  does  not  depend  on  t. 


11m 

h+0 


P0(t+h)  -  P0(t) 


P  (t)[-X  h+o,  0(h)) 

11m  -2 - 2 - - 

h-K)  h 


11a  PlCt+h)-p1(t)  .  lin  Fi(t)[-‘Xih4o2t1(h)]-tfi,1(t)[li,1hto1  iizj0Q]+o^ 
h-*0  h  h+0  h 

-  -AjP^t)  +  Xj^Pj.^t) 

so 

*  -*0po(t> 


and 


Pj'(t)  -  -XjPjO:)  +  Aj„1Pj_1(t). 

Solve  this  infinite  set  of  differential  equations  with  the  initial  con¬ 
ditions  by  using  the  integrating  factor  e^j*  for  the  set  _>  0 

PQ(t)  -  e-*0t, 

p3(t)  -  Xj-i  .'XJC  «XJ*  Pj-1  (Jt)  dx,  j  -  1,2 . 

Stating  the  postulates  of  a  pure  birth  process  in  terms  of  aircraft 
accidents  produces  the  following  assumptions: 


(1)  Of  the  Q|(  accidents  occurring  up  to  accumulated  time  tfe, 


the  last  one  occurred  exactly  at  time  t^. 

(2)  The  occurrence  of  accidents  follows  the  pure  birth  pro- 
cess  with  rate  An,  n  *  0,1,2,  ....  If  Pn(t)  denotes  the  probability 
of  exactly  n  accidents,  n  *  0,1,2,...,  by  accumulated  time  t,  then 

<*>  %'w  •  ar  po(t)  *  “W0 

(11)  P„'(t)  -  3F-  pn<c)  '  -1npn<c)  +  VlVlW' 

where 

P0(O)  -  1 
and 

Pn(0)  -  0,  n  >  1. 


After  examination  of  several  plots  of  n^  versus  tfc  e(tk))  for  air¬ 
craft  accident  data,  the  functional  relationship 

An  *  a  +  yyn,  n=0,l,2,... 

was  chosen  to  express  the  relationship  between  n^  and  A^* 

If  a  >.  0,  y  >  0,  and  0  <  y  <  1,  then  u  is  the  limiting  rate;  n  +  y 
is  the  initial  rate;  and  yn  is  the  fractional  part  of  y  remaining  after 
the  nth  accident. 

If  p  *  1,  then 


An  =»  a  +  y 

*  X, 

A  being  a  constant.  This  is  the  case  where  the  number  of  accidents  up 
to  accumulated  flight  tisc  t  has  the  Poisson  distribution  with  parameter 
At,  and  the  accumulated  flight  time  up  to  accident  number  n  has  the  gamma 
distribution. 

If  y  >  0  and  y  >  1,  the  number  of  accidents  as  a  function  of  t  will 
blow-up,  i.e.  there  will  be  a  positive  probability  that  the  number  of 
accidents  will  become  infinite  in  finite  time. 

If  a  -  0,  y  >  0,  and  y  <1,  the  rate 

An  *  YPn 

will  approach  zero  as  a  limit,  but  the  number  of  accidents  as  a  function 
of  t  will  be  unbounded  (logarithmic  increase). 
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II.  ESTIMATION  OF  THE  ACCIDENT  RATE  PARAMETERS 


For  data  such  as  chat  indicated  in  Table  2,  when  the  accident  rate 
between  the  nth  and  n  +  1st  accidents  is  taken  to  be  of  the  form 


Xn  =  a  +  ^n’  n  “  °»1»2 . 

the  computations  described  in  this  section  will  result  in  a  least  squares 
estimation  of  the  unknown  parameters  a,  y ,  and  p.  The  choice  of  a  quan¬ 
tity  Q  =  Q(a,y,p)  to  be  reduced  to  a  minimum  by  the  solution  a,  y,  and  p 
is  motivated  by  the  formula 

nk~l  , 
c(Tk)  -  Z  —j— 

J"nk-i  i 

discussed  in  section  I,  and  a  desire  for  a  close  fitting  cumulative  num¬ 
ber  of  accidents  versus  cumulative  fright  hours  curve  to  the  data  points 
of  Figure  2. 


The  least  square  estimates 

which  minimize  the  quantity 

m  nfc-1 
Q  *  h  l  (  l 
k-1  j-0 


m 


of  a,  y,  and  p  will  be  taken  to  be  those 


The  minimizing  a,  y,  and  p  are  to  be  found  by  refining  some  initial 
estimates  aQ,  yg,  and  pQ  by  Newton-Raphson  iterations.  Initial  estimates 
may  be  obtained  by  examining  a  nk  versus  observed  accident  rate  Nk/Tk 


17 


(see  Table  2)  scatter  diagram.  With  approximate  smooth  values  of  the 
accident  rate  for  »  0,  ■  some  Intermediate  value,  and  nfc  ■  ® 

(assuming  y  <  1),  solve  for  aQ,  y0,  and  yQ. 

The  minimising  a,  y,  and  y for  Q  satisfy  the  equations 


3  a  3  y 


3  y 


0. 


On  the  other  hand,  for  points  (a,  y,  y)  generally,  in  a  neighborhood  of 
(  a0,  yQ,  yQ),  the  Taylor  expansion  of  about  the  point  (aQ,  yQ,  yQ) 
is  given  by 


+  (y  -  y0)  C3a3y)o  +  higher  order  terms, 

where  (  )n  denotes  the  value  of  the  function  enclosed  within  the 

3Q  30 

parentheses  at  the  point  (aQ,  y0,  y0).  gy  and  -yj*  have  similar  expan¬ 
sions.  (Formulas  for  all  of  the  first  and  second  order  partial  derivatives 
of  Q  are  found  in  an  appendix).  Setting  the  left  side  of  each  of  the 
expansions  to  zero  (since  a  point  with  being  sought), 

dropping  higher  order  terms,  and  rearranging  what  remain,  lead  to  the 
system  of  linear  equations 
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(a  -  aj< 


3a3y 


>0  + 


(y  -  y„)C  7"?  +  ~4  >f 


3y3y 


3p" 


-<*>. 

9u  1 


whose  solution  (aj,  y^,  Wj)  should  be  an  inproved  estimate  of  the  point 
where  each  of  the  first  order  partial  derivatives  of  Q  vanishes. 


When  the  above  computations  are  repeated  with  a^,  y^>  and  y  j  re¬ 
placing  a0,  Y0f  and  yQ  wherever  they  occur  in  the  description,  a  third 
approximation  (a2,  y2,  P2)  is  obtained.  The  process  can  be  continued 
indefinitely.  The  sequence  of  approximations  (a^,  yi»  Pi)  for  i-0,1,2,..., 
converges  to  some  (a,  y,  IT)  which  satisfies  the  equations 


32  -  32  -  32  -  o, 

3a  3y  3p 

provided  the  initial  approximation  (a0,  y0,y0)  is  sufficiently  close  to 
(a,  Y,  p). 


In  the  computer  program  written  for  the  above  computations,  the 
system  of  linear  equations  is  considered  to  be  in  the  variables  (a  -  a^), 
(y  -  y^),  and  (p  -  p^)  (for  the  ith  iteration  where  (  )i  replaces  (  )0 
for  the  coefficients  of  the  equations)  and  the  system  is  solved  by 
Cramer's  method.  Then  a^,  y^  ,  and  y^  are  added  to  the  solutions  to 
obtain  a^+p  y^+p  and  y^+^.  The  partial  derivatives  of  Q  vanish  at 
externals  other  than  (local)  minima.  To  insure  a  correct  solution  set, 
the  program  prints  out  a^,  y^,  y^,  and  after  each  iteration  for  the 
programmer's  inspection.  If  c»i+i»  Yi+i,  Pi+i  yield  a  Q^+i  <.  Qj,  the 
program  continues  to  iteration  i+2.  If  >  Qj,  or  Qi+i  “  Qi  ■  Qi-1, 

then  a  ■  aj,  y  ■  yi,  y  ■  y^  are  taken  as'  the  solutions.  However,  if 
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Qj  >  Q0,  i.e.  an  increase  in  Q  results  on  the  first  refinement  of  (aQ, 

Y0,  Uq) ,  an  error  message  is  printed  out  and  computations  are  halted. 

In  this  case  new  trial  values  for  ao,  yo*  mo  must  be  entered.  Solutions 
are  obtained  in  about  6  iterations. 

After  a,  y»  and  M  have  been  found,  the  computer  program  proceeds 
to  compute  and  print  out  the  expected  cumulative  flight  hours  to  accident 
number  n, 

n-1  i 

e(tn)  -  Z  — 
j«°  Xj 

and  the  accident  rate 


Xn  *  a  +  YUn 

for  n  *  n0,  2nQ,  3nQt  .  (nQ  is  an  input  to  the  program  which  is  chosen 

to  be  10  in  the  sample  case.)  Plots  of  n  versus  e(tn)  and  Xn  versus  e(tn) 
should  be  made  and  the  data  points  of  Table  2  should  be  superimposed  on 
the  plots  as  a  check. 
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III.  APPLICATIONS  AND  EXTENSIONS  OF  THE  THEORY 


The  method  of  analysis  that  has  been  presented  in  this  report 
could  be  used  to  evaluate  a  wide  variety  of  data  that  relates,  or 
appears  to  relate,  to  aircraft  accidents.  In  each  case,  cumulative 
flight  hours  (time)  would  be  used  as  the  independent  variable.  Various 
sets  of  graphs  can  be  compared  for  indications  of  trend. 

A  point  of  departure  for  a  first  analysis  of  a  set  of  accident  data 
can  be  the  cumulative  number  of  aircraft  damaged.  The  graphs  of  this 
preliminary  phase  for  several  different  models  of  aircraft  may  indicate 
other  types  of  evaluation  and  comparison  that  should  be  conducted.  For 
example,  various  models  of  fighter  (or  bomber  or  cargo)  type  aircraft 
may  be  compared  to  see  which  models  have  lower  accident  rates.  The 
conclusions  that  are  drawn  from  this  examination  may  lead  to  the  scruti- 
nization  of  still  more  specific  aircraft  characteristics,  such  as  single 
engine  versus  twin  engine,  afterburner  or  no  afterburner,  and  aircraft 
weight. 

The  initial  sets  of  data  can  be  broken  down  into  the  five  aircraft 
damage  classifications  as  specified  in  the  "Navy  Aircraft  Accident, 
Incident,  and  Ground  Accident  Reporting  Procedure"  (OPNAV  Inst.  3750. 6F 
of  15  March  1967) : 

(1)  ALFA,  alfa  damage  (destruction  or  loss) 

(2)  CHARLIE,  substantial  damage 

(3)  DELTA,  minor  damage 

(4)  ECHO,  limited  damage 

(5)  FOXTROT,  no  damage. 
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These  cases  may  point  out  something  about  the  nature  of  accidents  in 
which  a  particular  model  becomes  involved.  For  instance,  it  may  be 
indicated  that  an  aircraft  is  highly  prone  to  accidents  resulting  in 
limited  damage  (ECHO).  This  result,  in  turn,  may  warrant  the  inspection 
of  the  phase  of  operation  in  which  the  accidents  happen.  Phases  of 
operation  could  be  classified  in  the  following  manner: 

(1)  Static  (engines  running),  incident  to  flight 

(2)  Taxiing,  incident  to  flight 

(3)  Takeoff 

(4)  Inflight 

(5)  Landing 

(6)  Waveoff  (go-round) 

(7)  Nonflight. 

For  use  by  the  Navy,  it  may  be  of  particular  value  to  look  at  the 
statistics  relating  to  embarked  and  disembarked  aircraft.  These 
statistics  could  then  be  broken  down  into  questions  about  such  matters 
as  the  type  and  size  o 2  the  carrier  and  the  length  of  time  that  the 
carrier  has  been  deployed.  An  increase  in  accidents  as  the  time  deployed 
increases  may  indicate  an  increase  in  the  operating  rate  and  commitments, 
over-confidence  on  the  part  of  pilots  and  crew,  fatigue,  and  anticipa¬ 
tion  of  the  approaching  rest  and  recreation,  all  of  which  could  be 
termed  unmeasurable  factors.  Carrier  based  causes  which  could  be  checked 
might  include  differences  in  the  methods  of  carrier  operations;  carrier 
operations  personnel,  such  as  lSO  (landing  signal  officer),  training 
experience  and  rating;  pitching  of  the  carrier  deck;  and  CVA. 
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The  number  of  years  an  aircraft  has  been  In  service  and  the  length 
of  time  since  major  overhaul  could  present  some  interesting  graphs. 

These  could  bring  into  play  statistics  about  material  failure,  malfunc¬ 
tion,  quality  control,  and  maintenance  procedures.  Also,  effects  of 
the  use  of  special  equipment,  machinery,  and  aircraft  support  equipment 
utilized  on  the  aircraft  and  the  ground  base  may  be  worth  noting. 

Other  flight  related  variables  that  could  be  considered  are  the 
time  of  day  or  night  of  the  flight;  season  of  the  year;  weather  condi¬ 
tions  including  wind,  sea  state,  cloud  coverage,  ceiling  visibility, 
temperature,  and  dew  point;  length  of  time  in  flight;  and  flight  altitude 
and  airspeed  at  the  time  of  the  accident. 

The  pilot  being  such  an  important  variable  in  flying  may  justify 
comparison  of  data  based  on  facts  like  total  pilot  time,  years  of 
military  flying,  pilot  rank,  pilot  time  in  the  specific  aircraft  types, 
pilot  time  in  the  last  three  months,  and  night  pilot  time  in  the  last 
three  months. 

Another  class  of  accidents  that  occurs  involving  aircraft  (not 
incident  to  flight)  is  ground  accidents  with  non-aircraft  vehicles. 

Areas  for  study  in  this  situation  might  include  causes  like  Improper 
action  by  the  operator,  mechanical  failure  of  equipment,  improper 
operator  action  accompanied  by  mechanical  failure,  unforeseen  occur¬ 
rence,  and  personnel  other  than  the  equipment  operator. 

By  plotting  several  related  curves  using  accumulated  aircraft  hours 
as  the  axis  of  abscissa  will,  hopefully,  reveal  some  relationships 


between  people  and  other  aircraft  accident  factors.  These  results  may 
then  be  used  in  decisions  about  matters  such  as  aircraft  selection,  size 
of  aircraft  inventory,  airframe  spares  procurement,  pilot  selection,  and 
carrier  operation  Methods . 


IV.  ALTERNATE  APPROACHED  TO  THE  PROBLEM 


The  analysis  of  aircraft  accident  information  is  of  great  and 
immediate  importance  to  the  U.  S.  Navy,  as  well  as  the  U.  S.  Air 
Force.  Studies  are  constantly  being  conducted  to  find  ways  tc  evaluate 
the  available  data  in  order  to  come  to  conclusions  and  decisions  con¬ 
cerning  non-combat  aircraft  accidents. 

Some  of  the  studies  and  evaluations  conducted  ar t  discussed  here 
to  show  other  approaches  that  have  been  used  in  the  analysis  of  aircraft 
accident  data. 

The  use  of  plots  of  cumulative  accidents  as  they  occurred  in  time 
in  aircraft  accident  analysis  is  well  known.  If  the  curve  for  the  given 
values  closely  approximates  a  straight  line,  the  accident  rate  is  con¬ 
stant  (see  Figure  3).  When  the  accident  rate  is  a  constant,  i.e.,  the 
hazard  function  given  by: 


h(t) 


f(t) 

1  -  F(t) 


is  constant,  say  h(t)  -  X,  then 


f(t)  -  h(t)  [1  -  F(t) ] 

■  h(t)exp[  -  jjj  h(t)  dt]  . 


Denote  the  slope  of  the  graph  by  X,  so  the  density  function  becomes 


f(t)  represents  the  exponential  probability  function  where  A  is  the  rate 
of  accident  occurrence . and  1/A  is  the  mean  time  to  occurrence. 

The  maximum  likelihood  estimator  for  the  sample  occurrence  rate 
is  normally  distributed  about  the  true  rate,  since  the  distribution  of 
the  sample  mean  is  asymptotically  normal  as  the  sample  size  approaches 
infinity.  A  confidence  interval  about  the  true  mean  is  needed  since 
the  true  variance  is  not  known.  Using 

t  -  ~  *) 

S 

A 

where  t  is  distributed  with  (n-1)  degrees  of  freedom,  A  is  the  sample 
mean,  S  is  the  sample  standard  deviation,  A  is  the  true  accident  rate, 
and  n  is  the  number  of  accidents }  then 

A  *  t  a/2  •  S 

StPT 

gives  a  (l-o)  confidence  interval  for  the  true  aircraft  accident  rate  A. 

Another  method  of  attacking  the  accident  prediction  problem  in¬ 
volves  the  systematic  invest igation  and  evaluation  of  statistical  infor¬ 
mation  in  order  to  derive  a  t' ingle  equation  (or  perhaps  a  set  of  equations) 
which  can  be  used  to  estimate  aircraft  attrition. 

This  is  the  process  used  by  T.E.  Anger  in  The  Estimation  of 
Peacetime  Aircraft  Attrition,  a  Center  for  Naval  Analyses  research 
contribution.  The  variables  Mr.  Anger  felt  affected  total  attrition  of 
forces  of  aircraft  included; 

(1)  Total  flying  hours  of  aircraft  force 
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(2)  Proportion  of  total  flying  done  from  carriers 

(3)  Empty  weight 

(4)  Maximum  speed 

(5)  Number  of  engines. 

From  the  results  of  statistical  measures,  the  variables  listed  above 
are  incorporated  into  an  attrition-estimating  equation  (the  actual 
equation  is  classified). 

W. E,  Mooz  assumes  that  the  cumulative  number  of  aircraft  destroyed 
is  a  function  of  the  cumulative  number  of  flying  hours  (dependent  varia¬ 
ble)  in  his  Rand  Corporation  Memorandum  Relationships  for  Estimating 
Peacetime  Aircraft  Attrition.  This  function  appears  as  a  reasonably 
straight  line  on  log-log  paper,  and  therefore  demonstrates  that  the 
cumulative  number  of  aircraft  destroyed  is  a  decreasing  exponential  func¬ 
tion  of  the  flying  activity  of  the  aircraft.  Mooz  states  this  is 
"evidence  that  the  attrition  pattern  of  a  given  type  and  model  of  air¬ 
craft  was  subject  to  an  orderly  learning  process  which  continues  through¬ 
out  its  flying  life".  The  exponential  function  y=ax^  is  used  to  approximate 
the  accident  curve. 

According  to  this  study  conducted  by  Rand,  the  quarterly  attrition 
rate  is  not  as  useful  a  variable  as  many  people  think.  Quarterly  attri¬ 
tion  rates  do  not  involve  equal  flying  times,  and  the  number  of  hours 
flown  varies  for  the  different  models.  Since  flying  programs  differ, 
there  can  be  no  true  comparison  between  quarterly  attrition  rates,  and 
also  attrition  does  vary  over  the  flying  life  of  the  aircraft.  To 
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sunmarize,  Mooz  says  that  "attrition  rate  is  awkward  £rom  a  statistical 
standpoint" . 

A  general  conclusion  is  made  by  A.  J.  Gross  and  Milton  Kamins  in 
Reliability  Assessment  in  the  Presence  of  Reliability  Growth.  They  show 
that  there  is  no  evidence  that  would  lead  to  the  generalization  that  the 
age  of  any  fighter  should  result  in  any  increase  in  accident  rates,  pro¬ 
vided  normal  preventive  mainf'  -.ace  and  product  improvement  are  continued. 

Three  statistical  examinations  of  accident  and  attrition  data  for 
jet  fighters  is  presented  and  compared  by  Milton  Kamins  in  Jet  Fighter 
Accident /Attrition  Rates  in  Peacetime:  An  Application  of  Reliability 
Growth  Modelling.  In  the  first  model,  the  reliability,  Rk,  at  any 
stage  k  of  the  process  is  given  by 

Rk  "  R.  '  f 

where  R^  is  the  ultimate  reliability  and  c  is  the  total  amount  of  relia¬ 
bility  growth  that  can  be  achieved  from  stage  1  to  stage  infinity  (i.e., 

-  Rj).  Accident  rate  (unreliability)  can  be  treated  as  the  complement 
of  accident  reliability  which  is  very  close  to  1.0.  Accident  rate,  F, 
is  expressed  as  events  per  hundred  thousand  flying  hours  or  hundred 
thousands  landings. 


i  -  Fk  -  1 


F 

00 


£ 

k 


where 


c  -  F  -  F 

1 
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so 


rk-F-+t 


Let  F  +  0 

CD 

log  *  log  c  -  log  k 
and 

c  -  F!  . 

Therefore 


log  Fk  -  log  Fj  -  log  k. 

Maximum  likelihood  estimates  are  then  developed  for  the  parameters  F^ 
and  c  in  this  hyperbolic  model  of  reliability  growth.  For  cases  eval¬ 
uated  by  this  model,  landings  rather  than  flying  hours  are  used  as  a 
measure  of  exposure,  since  they  appear  to  be  a  better  indicator  of  risk. 
Where  material  failure  is  a  contributing  cause  to  accidents,  the  accident 
rate  seems  to  be  a  function  of  the  number  of  years  in  service.  In  order 
to  cause  the  model  to  appear  as  a  straight  line,  the  independent  variable 
is  the  reciprocal  of  the  number  of  years  in  service  (see  Figures  4A  and 
4B).  With  this  type  of  graph,  comparisons  can  easily  be  made  between 
aircraft  which  involve  carry-over  technology  and  experience.  Comparison 
of  these  figures  A  and  B  show  that  even  though  the  ultimate  accident 
rate,  Fw  ,  is  near  zero  for  both  aircraft,  it  is  sufficiently  higher  in 
B  than  in  A  to  cause  the  curves  to  cross  at  about  the  8-year  point, 

Thus  it  is  shown  that  even  though  there  is  improvement  with  calendar 
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time,  the  second  aircraft  becomes  more  prone  to  material  failure  accidents. 
Other  comparisons  can  be  made  to  indicate  material  reliability,  such  as 
the  fraction  of  all  accidents  in  which  material  failure  is  a  factor  by 
years.  The  same  graphs  and  evaluations  can  also  be  made  for  all  causes 
of  accidents. 


The  second  model  of  reliability  growth  examined  is  the  (negative) 
exponential  model  in  which  the  reliability  Rk  at  stage  k  of  a  develop¬ 
ment  process  can  be  expressed  as: 

Rk  *  1  -  ae'0k 

a  represents  the  total  amount  of  growth  that  can  be  achieved  from  stage 
zero  to  stage  infinity,  and  B  is  a  measure  of  the  rate  of  growth. 

A  third  model,  called  a  learning  curve,  states  that  the  total 
number  of  accidents  (or  losses)  Lx  depends  on  the  total  usage  x  (flying 
hours  or  landings) . 


Ljj  -  AxB 

where  A  is  the  risk  for  the  first  flying  hour  or  landing,  and  B  is  the 
complement  of  the  rate  of  learning.  The  parameters  are  estimated  by  a 
weighted  least-squares  technique. 

Ramins’  three  models  are  compared  on  the  basis  of  a  chi-square  good¬ 
ness  of  fit  test  of  the  observed  data  by  year  against  the  predictions 
of  each  model  with  parameters  estimated  from  the  data.  Results  of  the 
evaluations  and  comparisons  of  the  three  models  show  that  the  hyperbolic 


33 


model  of  reliability  growth  is  preferable  for  three  reasons.  It  seems 
to  represent  trends  more  accurately,  is  easier  to  use,  and  useful  con¬ 
fidence  limits  for  past  experience  or  future  projections  of  accident  or 
attrition  rates  can  be  calculated. 

John  M.  Cozzolino  deals  with  the  infant  mortality  effect  of  statis¬ 
tical  reliability  theory  in  a  paper  titled  Probabilistic  Models  of 
Decreasing  Failure  Rate  Processes.  Infant  mortality  implies  a  decrease 
of  the  conditional  probability  of  failure  of  a  device  with  increases  in 
age.  Cozzolino  compares  his  four  models  of  decreasing  failure  rate  pro¬ 
cesses  based  upon  the  population  heterogeneity  hypothesis  and  incorporating 
explicit  repair  assumptions  (see  Table  3). 

The  various  models  presented  could  be  used  in  predicting  different 
aspects  of  aircraft  accidents.  For  instance,  the  component  variability 
model  could  be  used  in  predicting  the  aircraft  accidents  which  will  re¬ 
sult  in  total  destruction  (non-repairable). 

All  the  methods  and  models  that  have  been  described  have  some 
theoretical  application  to  the  case  of  aircraft  accident  prediction. 

Bach  has  its  advantages  and  disadvantages,  so  it  is  necessary  to  choose 
a  method  (or  methods)  that  most  accurately  fits  the  situation  being 
evaluated. 
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COMPARISONS  OF  FOUR  MODELS  OF  DECREASING  FAILURE  RATE  PROCESSES 
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V.  COMPUTER  PROGRAMS 


The  following  are  listings  of  the  two  programs  used  in  this 
analysis. 

The  first  program  finds  a  linear  solution  to  a  set  of  equations 
by  Cramer's  Rule.  The  second  program  modifies  the  input  data,  mini- 
mizes  Q,  and  determines  the  parameters  a,  y ,  and  p. 


36 


lidear  program  urTir:n 


sunnnuTinE  lipear 

COMIlOr  A(3^4)^S( 3) 
nif!Ef!5l 0!’  B(3,4),DET(4) 

DO  u  1=1,3 
DO  i:>  J=l,4 
W  DCI/J)  =  A(  I , J) 

DC  GO  C=l,4 
DO  2 0  1=1,3 
2C  CO,!.*)  =  A ( 1,4) 

I  =  1 
J  =  2 
K  =  3 
R  =  1. 

DET(fl)  =  0. 

DO  53  L=l, C 

DET(!I)  =  PET(f!)  +  R*B(l, I )*B(2,J)*C(3,K) 
JJ  =  J 

ID  (R)  3D, 33,M-' 

3D  J  =  K 
K  =  JJ 
no  TO  50 
4  C  «J  —  I 
!  =  JJ 
5D  R  =  -R 

DO  60  1=1,3 

en  c(i,i:)  =  AU,n) 

no  70  1=1,3 

70  SCI)  =  DE.T(  |  )/DET(4 ) 

RETURf: 

END 


MAIN  PROGRAM  LISTING 


i 


I 

I 

i 


j 

\ 

1 

( 

{ 


COIUIOM  QAA,DO.AG,  DRAM,  PAG,  P.GG,  DOOM,  PAM,  PGM,  PI  Hi,  SA,SG,  SM,  ALA, GAM,  MU 
DIMERS  I  OH  ACT(100),Cf!(100),CT(10('),E(100),EA(100),EAA(106), 
1EAG(100) ,  EAM(  10[!),EPi(100),EGG(lC0),EGIK100),  EGM1 (100 )  ,EGi!2  (100) , 

2 EM(  10  0 ) , EMM ( 1 0  0 ) , FMM1 ( 10  0 ) , EMM2 ( 1 0  0 ) ,  RATE  ( 1 0 0 ) , T ( 1  fJ 0  )  , T  I  ( 1  [* 0 ) 

REAL  I  ACT,  I  LAM,  JMJ,  JI1J2,  J1JMJ,  J2I1J2,  LAC,  LAJ,LAJ2,LAJ3, LAM( 100 ) , 
1LAO, LMC,MC,MiJ,f1J2,MO,MU,Ml, f!( 100), Ml  (100) 

DO  460  HCASE*1, 4 
WRITE  ( G  ,  1 0 )  RCASE 
10  FORMAT  (////, 32”, 5MCASE  ,11) 

TEII4  -  10. **(-4) 


READ  (5,20)  IM 
20  FORMAT  (1C) 

I  =  1 
TT  =  0. 

DO  90  <J=1,  IM 

READ  (5,30)  T(J),II(J) 

30  FORMAT  (2F10.0) 

T I ( J )  =  T(J) 

Ml ( J )  =  N(J) 

IF  (IK J) )  40,70,C0 
40  WRITE  (6,50) 

50  FORMAT  (/ , 2SI!  11(10  SHOULD  HOT  3E  NEGATIVE) 
WRITE  (6,60)  K,IKK) 

60  FORMAT  (10X,4IIK  =  ,  I  6, 4X,  7IIIK 10  =  ,F6.0) 

GO  TO  460 
70  TT  =  TT+T( J) 

GO  TO  90 

80  T( I )  =  (TT+T( J) )*TEN4 
NCI)  =  l!(  J) 

RATE(I)  =  fl(  I  )/T(  I ) 

I  =  1+1 


‘  TT  -  0. 

90  CONTINUE 
M  =  1-1 

;  CT( 1)  =  T( 1 ) 

C  N ( 1 )  =  N(l) 

DO  100  I  =  2  ,  M 

CT( I )  =  CT( I -1)+T( I ) 

100  CN(  I  )  =  Cf! (  I  -1)+N(  I  ) 

WRITE  (6,105) 

105  FORMAT  (//,10X,10IIINPUT  DATA,  17X,  25IIMOD I  F  I  ED  CUMULATIVE  INPUT, 
15H  DATA) 

WRITE  (6,110) 

110  FORMAT  (/,4X,  GHNUMDER,  8X,  8IIA I  RCRAFT,  8X,  6HNUMBER,  5X,  9IIA I  RCRAFT  , 
15IIH0URS,  4X,  8HACCI  DENT) 

WRITE  (6,115) 

L15  FORMAT  (1311  OF  ACCIDENTS, 7X, 5IIII0URS,  6X,  12II0F  ACCI DENTS, 4X, 

110HX  10**(-4),8X,  4IIRATE) 

DO  120  1=1, M 

120  WRITE  (6,130)  fl I  ( I ) , T I  ( I ) , CM ( I ) , CT(  I ) , RATE ( I ) 


1^3i,ii1sJi+1 

DO  135  !»::r:i,|M 
135  WRITE  CG,13')  r:i(|),TI(|) 

137  IF  (M-3)  140,  lGp, 160 
140  WRITE  (6,15i.)  i: 

!50  FORMAT  (/,3',:i:  MODIFIED  IIIPUT  DATA  WAS  Of’LY  II 
12911  PEP.ICDC  Will  Cl!  IG  LEST  TI'Af*  *)  #  ' 

GO  TO  4  GO 
ICt  WRITE  (6,170) 

17t  FORWAT  (/,  52!!  IMPUT  INITIAL  AMD  FI!!AL  RATfG  AMD  !  MTEPMFD  I ATF 
1151!  At!!)  CUMULATIVE)  A,u  *•  l  ti.ilhu I ATL 

WR I  TP  ( G, 175 ) 

175  FORMAT  (29M  ACCI  DFf'TS  USIUC  FORMAT  4F5  1) 

REAP  (9, 1C  ')  LAO,ALO,LAC,Cf!A 
ICC  FORMAT  (4F5.1) 

GAO  =  LAC-ALP 
MC  =  ( LAO.-AL°)/PAO 
EMC  =  ALOC(MC)/CUA 
MO  =  EXP(L,MC) 

write  (0,2:0) 

m  ro^ru^^ncc^vr  approx  i  i:ati  or:E ) 

22'i  KXI  J/'ri;'5,'AL,''l!A,U'X,5l'nAMr:A,li:<,S,:;iM(12X,lHO) 

JPASS  ='  o 
AL1  =  o. 

CA1  =  0. 

Ml  *  0. 

Q  =  932959999.  4 

230  Q1  -  Q 

DO  2G0  K=1,M 

J  =  -1 

E(K)  =  0.  ^Qh> 

EA(  K)  =  0. 

EG(K)  =  0.  ^/k 

EM(K)  =0. 

EAA(K)  =  0. 

ECGOO  =  0. 

EMMl(K)  =  0, 

EMM  2  (10  =  0, 

EAO(K)  =  0. 

EAIKIO  =  0. 

EGIIl(K)  =  0. 

EGM200  =  0. 

240  J  =  J+l 
MJ  ®  M0**J 
MJ2  =  MJ*MJ 
JMJ  3  «J*MJ 
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!  5  [> 


Jf'J2  =  J*MJ2 

J2CJ2  =  j*ji:j2 

JlJfJ  =  <  J-l )  *  Jf  M 
LAJ  =rfALn+nAO.*MJ 
LAJ  2^  =  LAJ*  LAJ 
LAJ 3  =  LAJ*  LAJ 2 
fl(  K>  =  !:( K)  +  l .  /  LAJ 
HA(K)  =  FAd>)  +  l  ./LAJ  2 
EG(K)  ■  C0(K)+.*:J/LAJ2 
nr(io  =  Ef:(io+j:u/LAj:.' 

rAA(r)  =  FA  A  (!')  +  !./  LA  J  :> 
EPP(K)  =  EGC(lO+rj2/LAJ3 
CMI’KIO  =  Fi;;il(r,)+J2.MJ2/LAJ3 
ni'.r,2(K>  =  C!:M2(K)  +  J1J!M/LAJ2 
FAPd'J  =  FA0(I/J+:,J/LAJ3 
CAf  1(10  =  FAr(K)+JflJ/LAj:'» 
EdlldO  =  Fd.l(K)+J.,:J2/LAJ3 
EfSf'2 ( r:)  EC."2(  K)+Jf'J/LAJ2 

IF  (J-ClKl'J  +  l)  241/  2F|-/  240 
E ( K )  =  r(K)-CT(K) 

FACO  =  -EACO 

rrc i )  =  -eg( k) 
f;i(  i')  =  ( -nAo/'u' ; 

FAA(K) 

EGr(K) 


NCT  REPRODUCIBLE 


2  Cl' 


EfV'dO 
FA^(r) 
FAMd'.) 
FPIldO 
0  =  p. 
PA  =  f!. 
OP  =  <1. 

on  =  r' . 


2  .  *  FA  A  ( T. ) 

2 . *  FPP( K ) 

( (2.*PAP*GAP)/(f!0*n0))*FiridO-(PA0/(nC*l'n)  )*E!1M2(K) 
2 .  *FAr(  10 

( (2.*P.AP)/nn)*rAM(K) 

( ( 2 .  *PAC )  /no )  *  EGM1  (K)-(l.  /IIP )  *  cn  !2 1 


QAA 

npp 

Pill! 

PAP 

PAM 

PPM 


0. 

0. 

0. 

0. 

0. 


2/ 


no  2.::  ii 

p.  =-  n+Ed:)*r( id 

qa  =  PA+F(i:)*FA(P) 

rn  =  pn+rdd*FPd;) 

pi:  =  pii+f(i:)*!;i:(i:) 

paa  -  r  aa+ r  a  do  *  r  a  (  k  )  +  r  (K)  *  FA  a  (  k  ) 

C.GP  -  0PP+CP(i;)*FP(iO  +  F(K)*FGn(K) 

niii-  =  pnn+r;i:(K)*rM(i;)+r(id*FMn(id 
pap  =  PAp+EACid  *rPd;-)  +  F(K)*FAn(K) 
PA,";  =  PAI!+EA(i::)*n,dC)  +  r(K)*EAi;(K) 
ppi;  -  ppi;+FP(K)*n:(iO  +  F(K)*EPf:(id 
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0 


P  =  P.S*C 

i  ff'n  =  irrn+i 

i;piTr  (c#rco  At.r,r/\o,rr,r 

2Cf  fprpat  ( n r i r . r ) 

ir  (r-ri)  3;{,#2nr,3r:n 
29f  JPACS  =  JPASr+l 

ir  (A^T  ( A  L1-AI.0 )  + Arf/(PA1-PAP  )+ATC  (f’l-f'P ) )  3P0,31P,3&3 
sr.f  ir  (jPArr.-3)  3»:r^ 31j:,3m: 


33  0 

IF  (1 

rpp-3 )  320 

32r 

A I  P  = 

ALP+TPPL 

p«r  = 

PA:'+Tr‘!i: 

pp  = 

PP-T  -’’L 

pp  TP 

n  *  ft 

33;.: 

J  PASS 

=  0 

3  hi 

FA  = 

-pa 

sr  = 

-rr 

SP  = 

-rp 

DCAP 

=  pap 

PfAP 

=  PAP 

PPPP 

=  ppi: 

PALL 

LIPPAP 

A  LI  = 

ALC 

PA1  = 

GA  C 

P.l  = 

nr 

ALC  = 

*  LP+ALA 

PAP  = 

PAP + PAP 

P.P  = 

f’P+.PP 

rp  TP 

2  3i; 

3Pj> 

IT  (lfTT-2)  3 or 

3rr 

ALA  = 

!  A  LI 

PM*  = 

:  P'U 

M  j 1  = 

Ml 

1  L  V' 

=  a  l  a+p  v 

1  APT 

=  1 .  /  l  L  A  P 

r:p  = 

pp(p)  +  ir'. 

vA> 


K  ~  1 

i  =  ir 

nC  3  ."ft 

I  LAP  =  A  LA  +  n  AP*PP** ( J -1 ) 

I /PT  =  I  ArT+l ./ I  L.V 
IF  (J-L)  3  F  f  /  3  7  P  /  3  ■"  f‘ 

37P  I.AP( !". )  -  !  LAP 
APT (  Y. )  =  I  APT 
V  =  K  +  l 
L  =  L+lft 
3CU  PlTTIPL'P 

i.riTr  (c,3'ji:) 

3 ‘J (.  r PPMAT  ( / /  ,  1 9/':,  1 1! ’OUTPUT  DATA,/) 
l.TITF  (G,liftO)  A I  A  /  PAP,  PU 
4CC  F  PITA  T  (OP  ALPI'A  =  ,  IT  2 . 1, 2X,  EITAfTA 


=  ,iT2.c,2:;,5!ru  =  /ri2.o 


/DATA 


WRITP  (G,»il(3) 

FORMAT  (/,UX,10[T!JMULATIVr) 
l.'RITP  (C,412) 

FORMAT  (4X,Cmif!rFn,rX,UMAIRrRArT  FOURS ,  bX,  CKArr  f  nrr.j) 

WRIT*  (6,415)  ' 

FORMAT^nr  or  ArrmFrTS,7x,i0i»x  ic**(-4>,sx,4i\TATr) 

L  ■  10 

DC  430  d«l,K 

WRITF  (G,42fJ)  L,AOT(«J),  LAM(J) 

FCPI1AT  (2X,IG,12X,r9.4,5X,F9.4) 

L  *  L+lp 
no  TO  4G0 
WP.ITr  <G,45C> 

ccmwuF/#4r,!  TnY  Ari0TI’FR  MFT,,°D  ron  oPTAinr.T  intial  valutd 

STOP 

FOR 


tiOt 


VI.  SAMPLE  TEST  CASES 


Since  Navy  accident  data  is  classified,  it  has  been  necessary  to 
set-up  sample  statistical  information  that  will  serve  as  a  set  of  rep¬ 
resentative  cases. 

Included  in  this  section  are  a  listing  of  the  accident  data  used, 
the  computer  print  out  for  the  four  sample  cases,  and  the  graphs  that 
relate  to  each  case. 
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e 


INPUT 

DATA 

MOPIFIFP 

CUMULATIVE  INPUT 

DATA 

nuMP^n 

AIRCRAFT 

f'UMRFR 

AIRCRAFT  HOURS 

ACCIRFNT 

or  ACC  I  PELTS 

HOURS 

of  accidhits 

X  1(5**  (~4  ) 

RATF 

0. 

0. 

4. 

0.1C93 

23.6267 

If. 

1093 . 

9. 

0 .4230 

19.C4C4 

5. 

2545. 

12. 

0.774C 

8.5470 

3. 

3510. 

15. 

1.2109 

6. 8792 

3. 

4301. 

19. 

1.7314 

7.6849 

I*. 

5205. 

24. 

2C 3804 

7.7042 

5. 

0490. 

30. 

3.1803 

7.4451 

0. 

2059. 

35. 

4.2169 

4.8515 

5. 

IP 30C. 

39. 

4.9549 

5.4201 

if. 

7380. 

44. 

0.3804 

3.5075 

5. 

14255. 

51. 

7.7CC7 

4.977C 

7. 

14003. 

SC. 

9.2757 

4.7011 

7. 

14890. 

03. 

10.9013 

3.0758 

5. 

10250. 

or. 

12.5413 

3.048C 

5. 

10400. 

74. 

14.2587 

3.4937 

C. 

17174. 

79. 

10.2255 

2.5422 

5 . 

19CCC. 

85. 

18.2934 

2.9015 

6. 

2  0  0  7  G . 

91. 

20.1020 

3.2110 

C. 

Irene, 

99. 

21.8930 

4.6200 

O 

o  • 

17310. 

105. 

23.8840 

5.013C 

c. 

19910. 

115. 

25.9828 

4.7658 

10. 

20983. 

127. 

29. 9720 

3.0077 

12. 

39898. 

139. 

34.1742 

2.85C1 

12. 

42016. 

148. 

38.4804 

2.0871 

S. 

43122. 

163. 

43.4228 

3.0387 

15. 

49304. 

174. 

47.0901 

2.5777 

11. 

42073. 

185. 

52.0078 

2.5477 

11. 

43177. 

191. 

55.0704 

1.6355 

c. 

3G08C . 

201. 

209. 

59.9575 

2.3468 

10. 

42011. 

03.8697 

2.0345 

8. 

39322. 

223. 

67.4432 

3.9177 

14. 

35735. 

INPUT  INITIAL  AND  FINAL  RATES  AND  IKTERfJFPIATF  RATF  AND  CUMULATIVE 
-  accidfnts  usirm  format  4F5.1 
M.  0077  Ff’TFR  DATA. 

11.5  2.0  3.0  100.0 


k. 


^OI-REPKODUCI^^ 


CSSI VF  APPROX  I  TATI  OfS 


ALPHA 
2.00000000 
2 .13074589 
2 . 17228889 
1 . 05422230 


OAMMA 

9.5 0000000 
7.03795523 
0. 20930020 
4.97937107 


IV 

C. 97773002 
0 , 90045221} 
0.90200403 
0. 2CCC3957 


5C. 57025140 
10 . 520)31220 
10.00150752 
27.12139093 


OUTPUT  PATA 

ALPHA  *  2.17228009  rVUTA  ■  G.2093G02C  IMI  =  0.90208403 


rurtuLATivr 


fJUMRFR 

or  Acrmrrn 


A  I  RPRAFT  I'OURO 
X  l8**(-4) 


ACnPPHT 

RATF 


10 

1.2581 

7.5003 

20 

O  r  n  o  r* 

«.  •  v  o !  • 

0.0191 

30 

4.3010 

5.0,837 

40 

0.1(590 

5.2G99 

50 

0.1191 

4.7570 

00 

10.3350 

4.3301 

70 

12.7507 

3.9732 

00 

15.3881 

3.0754 

90 

10.2181 

3.4208 

180 

21.2405 

3.2193 

110 

24.4448 

3.04C2 

120 

27.8180 

2.9010 

130 

31.3482 

2.70,10 

148 

35.0194 

2.0804 

150 

38.0177 

2.5903 

100 

42.7291 

2.5202 

170 

40.7402 

2.4077 

100 

50.0385 

2.4100 

190 

55.0125 

2.3781 

200 

59.2519 

2.3440 

210 

03.5475 

2.3150 

220 

07.8911 

2.2919 

230 

72.2755 

2.2721 

ACCIDENT  RATE 


i 


FIGURE  8 

CASE  1 


APPROXIMATE  CUMULATIVE  ACCIDENTS  VERSUS  ACCIDENT  RATE 


0  80  160  240 

CUMULATIVE  ACCIDENTS 
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CAST  2 


IffPUT 

DATA 

MODIFIED 

CUMULATIVE  INPUT 

DATA 

NUfTFR 

AIRCRAFT 

NUMPFP 

AIRCRAFT  HOURS 

ACCIDENT 

OP  ACCIDENTS 

MOURS 

OF  ACCIDENTS 

X  10**(-4) 

RATE 

0. 

0. 

2. 

0 . 1693 

11.8134 

2. 

1093. 

5. 

0.4238 

11.7878 

3. 

2545. 

6. 

0*.  7748 

2.8490 

1. 

3510. 

7. 

1.2109 

2.2931 

1. 

4301. 

9. 

1.7314 

3.8425 

2. 

5205. 

12. 

2.3804 

4.0225 

3. 

0490. 

15. 

3.1863 

3.7225 

3. 

8059. 

17. 

4. 2169 

1.9406 

2. 

1030G. 

20. 

4.9549 

4.0050 

3. 

7380. 

22. 

6.3804 

1.4030 

2. 

14255. 

25. 

7.7CT7 

2.1333 

3. 

14003. 

29. 

9.2757 

2.6864 

4. 

14890. 

32. 

10.9013 

1.8455 

3. 

16250. 

34. 

12.5413 

1.2195 

2. 

16400. 

38. 

14.2587 

2.3291 

4. 

17174. 

42. 

16,2255 

2.0338 

4. 

19GG8 . 

45. 

18.2934 

1.4507 

3. 

20679. 

49. 

20.1620 

2.1406 

4. 

18686. 

55. 

21.8936 

3.4650 

6. 

17316. 

58. 

23.8846 

1 . 5008 

3. 

19910. 

63. 

25.CC2E 

2.3829 

5. 

20983. 

71. 

29.9726 

2.0051 

8. 

39898. 

78. 

34.1742 

1 .6000 

7. 

4201G. 

86. 

38.4804 

1.8552 

c. 

43122. 

94. 

43.4228 

1 . 0206 

8. 

49364. 

104. 

47.6901 

2.3434 

10. 

42673. 

111. 

52.0078 

1.6212 

7. 

43177. 

116. 

55.6764 

1.3629 

5. 

36686. 

123. 

59.9375 

1.0428 

7. 

42011. 

129. 

C3.8G97 

1.5259 

6. 

39322. 

138. 

67.4432 

2.5185 

9. 

35735. 

INPUT  INITIAL  AND  PINAL  RATES  AND  I NTERMFDI ATE  RATP  AND  CUMULATIVE 
*  ACCIDENTS  USINC  FORMAT  4E5.1 
M. 0077  ENTER  DATA. 

4.5  1.5  2.0  50.0 
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SUCCFSS1VF  APPROXIMATES 


ALPHA 

GAMMA 

MU 

0 

1.50000000 

3.00000000 

0.96479928 

34.07696533 

1.58477497 

2.00770378 

0.97043335 

15,92052460 

1.59876156 

1.94672012 

0.97252291 

11.428367G1 

3.484C4584 

4.59255791 

0.83409184 

249C , 21679GS7 

OUTPUT  DATA 

ALPHA  -  1.59876156  GAMMA  - 

1,94672012  MU 

«  0.97252291 

CUMULATIVE 

MUMPER 

AIRCRAFT  HOURS 

ACCIPERT 

OP  ACC  1  HERTS 

X  10**(-4) 

RATE 

10 

3. CISC 

3.1137 

20 

6.4653 

2.7453 

30 

10.3365 

2.4665 

40 

14.G02G 

2,2555 

50 

19.225G 

2,0958 

or 

24.1614 

1.9749 

70 

29.3G3S 

1 , BC  35 

80 

34.7879 

1.8142 

90 

40.3929 

1.7618 

100 

46.1432 

1.7222 

110 

52.0084 

1.6922 

120 

57.9638 

1.6695 

130 

63.9893 

1.6523 

140 

70.0690 

1.6393 
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CUMULATIVE  ACCIDENTS 


FIGURE  11 


CASE  2 

APPROXIMATE  AIRCRAFT  ACCIDENT  CURVE 


CUMULATIVE  AIRCRAFT  HOURS  X  10 
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ACCIDENT  RATE 


i 


FIGURE  12 
CASE  2 

APPROXIMATE  CUMULATIVE  ACCIDENTS  VERSUS  ACCIDENT  RATE 
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I 

I 


It! PUT  DATA 


i 


PUMPER 

OF  AOriDFfJTS 
1. 

0. 

2. 

7. 

4. 

4. 

1. 

1. 

3. 

2. 

3. 

4. 

4. 

4. 

4. 

4. 

5. 

5. 

5. 

7. 

5. 

13. 

11. 

12. 

6. 

15. 

12. 

13. 

3. 

7. 

11. 

0. 


A I Rf RAFT 
HOURS 
G4. 
3. 
130. 
334. 
G81 , 
1827. 
973. 
719. 

1430. 

1512. 

1726. 

2510. 

3290. 

3221. 

2247. 

2768. 

3385. 

4197. 

3990. 

4509. 

5512. 

9168. 

867G, 

10285. 

12499. 

11315, 

12271. 

15088. 

6869. 

3792. 

7083. 

5482. 


CAST  5 


MOMFirn  rUMULATIVF  irrUT  DAT A 


UUMPPR 
OF  ACrjpruTS 
1. 

3. 

10. 

14. 

18. 

19. 

20. 

23. 

25. 

28. 

32. 

36. 

40. 

44. 

48. 

53. 

58. 

63. 

70. 

75. 

88. 

99. 

111. 

117. 

13  2. 

144. 

157 ; 

1G0. 

1G7. 

178. 


aircraft  hours 
X  10**(-4) 
0.C0C4 
0.0197 

O. 0531 

0.1212 

P. 3P39 
0.4012 
0.4731 
0.61C1 

jj.  7673 
0.9399 
1.1909 
1.5199 
1.8420 
2. 0667 
2.3435 
2.F870 
3.1017 
3.5007 
3.9516 
4.5028 
5.4196 
6.2872 
7.3157 
3.5G5G 
9. 6971' 
10.9242 
12.4330 
13.1199 
13.4991 
14.2074 


ACriDFNT 

RATE 

156.2501 

150.3760 

209.5809 

58.7372 

21.8938 

10.2775 

13.9082 

20.9790 

13.2275 
17.3812 
15.9363 
12.1581 
12.4185 
17.8015 
14.4509 
14.7710 
11.9133 
12.5313 
15.5245 
9.0711 
14.1798 
12.6787 
11.6675 
4.8004 
13.2567 
9.7792 
8.6161 
4.3674 
18.4599 
15.5301 


ACCIOeM^oSAK1^1!8  AMD  RATE  AND  CUMULATIVE 

M.0077  FHTFR  DATA. 

42.0  10.0  12.5  60.0 
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success i  vf  approximations 

ALPHA  CAf'MA  Ml'  0 

10. C  C00C0CC  32. 00000000  0.95C3993C  0.59225416 

9 . 89008331  33.21765137  0.95824713  0.57824373 

9.86445332  32.5G52C184  0.05901471  0.57798862 

9.86572552  32.60165405  0.05897353  0.57799016 


OUTPUT  OATA 

ALPHA  *  9 .86445332  MAMMA  =  32.5652C184  f*U  =  0.95901471 

CUMULATIVE 


tmnrrp 

AIRCRAFT  IICUPF 

Acrpn- 

or  ArriprrTs 

X  10**(-4) 

nATF 

10 

0.2721 

32.2095 

20 

0.6347 

24.5684 

30 

1.0990 

19.5403 

40 

1 . 6685 

16.2315 

50 

2.3381 

14.0542 

60 

3.0954 

12.C215 

70 

3.9243 

11.6787 

80 

4.8082 

11.0583 

90 

5.7325 

10.6501 

100 

G.G856 

10.3814 

110 

7.G5S6 

10.2046 

120 

8.6451 

10.0883 

130 

9. 6407 

10.0118 

140 

10.G425 

9.9614 

150 

11.6483 

9.928,2 

160 

12.C5C9 

9.9064 

170 

13.GG72 

9.8921 

180 

14.6786 

9.C82C 

60 


FIGURE  14 


CASE  3 

ACTUAL  CUMULATIVE  ACCIDENTS  VERSUS  ACCIDENT  RATE 
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CUMULATIVE  ACCIDENTS 


TASr  4 


If! PUT  DATA 


MORI  FI m  rtJHJLATIVP  input  DATA 


or  AfrinFrTs  "KSF  or  S,:TS  "HSOT8  ArJ 

0*  c5*  b  0 • P531  37. 

8:  no*  f*  i4. 

2  *»•  0.31539  5. 

334.  r  fi  rin  ,• 

1  rn  i*  P.Cliil  3. 

,  Gol,  f,  (1  7(573  r 

1.  lf',7  7  W»/u/j  C, 

0.  973  1‘  0-93OO  5. 

1)  7io*  ®*  1.1909  7. 

X*  .?lb  10*  1.5199  3.1 

I.'  “•  l.S*2C  3.! 

i.  ijir  }*•  2-0OC7  c.i 

2  Hio*  }?•  2-3li5!:  J-' 

1.  3’or-  “•  2.0S2C  5.! 

1  H*  3.1017  2.: 

2  3.5007  5.1 

1  !;  •  ;•  3.0516  6.1 

2  24.  4.5028  3  r 

r  33"5’  29.  T.419C  it 

1.  1.107  *•  .-.tj-jo  3,»l 

„  34,  fi  0070  r  1 

2-  3990.  37  r*7 

3*  4509.  3r*  o’rjrc  2  •  J 

2.  5517  C.5C^6  0,5 

5  nirp*  42,  9.G971  3,5 

5*  orSb  45.  10.9242  2.4 

i  ic"?-  5  •  12.4330  3 ! 0 

1.  Izooo'  •  13-11M  2.9 

4.  1131*"  55*  13.4991  5.2 

3.  iilh:  n-Kn  7-° 

6*  15000. 

2.  Gr.r.9. 

2.  3792. 

7003 . 

0.  5402. 

.AmJrwsTiMMS'roorATl4rfiS  A,ID  l"TO,'iEn,/'TF  MTE  amd  cumulative 

H.C077  EWTCR  DATA. 

3.0  3.0  4.0  30.0 


mippp 

nr  ArnirrnTc. 


AIPFRAFT  f?C(JRS 
X  10** (-4) 
0.0531 
0.1212 
0.3039 
0.C161 
0.7G73 
0.9399 
1.1909 
1.5199 
1.2420 
2.0GG7 
2.3435 
2.G220 
3.1017 
3.5007 
3.9516 
4.5022 
5.419C 
6.2072 
7.3157 
2.5656 
9.G971 
10.9242 
12.4330 
13.1199 
13.4991 
14.2074 


ArriPFt.’T 
RATF 
37.GC42 
14.G243 
5.4735 
3.2031 
G.G138 
5.7937 
7.9681 
3.0395 
3.1046 
2.9002 
3.6127 
5.9084 
2.3827 
5.0125 
G. 6534 
3. 6284 
5.4538 
T . 7630 
2.9169 
0.8001 
3.5351 
2.4448 
3.9767 
2.9116 
5.2743 
7.0592 


•  ***■*.  jgi 


> 


sucpfssivf  approx imatiohs 


ALPHA 

3.00000000 

2.73181629 

2.71273613 

3.51560020 


GAMMA  MU 

5.00000000  0.94776571 

5.05465126  0.94653177 

4. 73335GG7  0.95164692 

5 . C42G530C  0.92134029 


OUTPUT  DATA 

ALPHA  -  2.71273613  GAMMA  *  4,78335667  MU  = 


CUMULATIVE 


N'UMPFR 

Al RrRAFT  HOURS 

ACCIOFHT 

OF  ACC| OFf’TS 

X  10**(-4) 

RATF 

10 

1.530G 

5.7748 

20 

3.5124 

4.5781 

30 

5.9293 

3.8491 

40 

8.7206 

3.4050 

50 

11.8037 

3.1345 

60 

15.0971 

2.9697 

70 

13.5333 

2.CG93 

Q 

3.32456112 
1 . C1308098 
1.53015137 
3.36095905 


0 . 95164G92 
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ACTUAL 


© 


© 

© 


CUMULATIVE 


FIGURE  18 


CASE  4 

CUMULATIVE  ACCIDENTS  VERSUS  ACCIDENT  RATE 


SUMMARY 


This  report  has  presented  a  method  of  predicting  aircraft  accidents 
using  the  pure  birth  process  and  relating  accident  rates  directly  to 
the  total  number  of  past  accidents. 

Other  methods  used  to  predict  aircraft  attrition  were  also  included 
in  order  to  demonstrate  weaknesses  and  strengths  in  these  various  methods 
and  in  the  pure  birth  method.  Then  depending  on  particular  characteristics 
of  a  case,  the  most  accurate  prediction  method  can  be  employed. 
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APPENDIX.  PARTIAL  DERIVATIVES  OF  Q 


Q  -  (<*»  Y»  y) 

m  nk-l  i  ® 

-  \  l  {  Z  7-  -  tk}2  -  ^  l  Fk2 
k-1  j-0  Aj  k-1 


nk-l  i 

where  Xj  *  a  +  yyJ,  Ek  *  l - tk. 

J-0  Xj 


92Ek  .  2E  HL 

9a9y  Xj3 


32£k  .li  j  j  yj 

9a9y  U  Xj3 

i!!k-  2x  J  -IeIhI 

9y9y  y  Xj3  m  Xj2 
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All  of  the  summations  indicated  above  run  from  j»0  to  j-nk-l.  The  sum¬ 
mations  below  run  from  k*l  to  k<*m. 


M  -  JE,  3Ek 
3a  1  Ek  -3oT 


3^  »  £  Ek  Bk 
3Y 


IS  .  r  „  i!k 

a..  ^ 


g.i[<£>*  +  *gi  i 


13  vr,  3Ek  ,2  32Ek  . 

3?"  Il<  WY  +EkW  1 

|f|.E1(iBL)a  +  Ek|^  , 


«(£><{£>♦*&) 


3a3y  3y3a 


"cFaSy 


ii£L.ii2  .  I((i*k)(i5S.)  +  Ek^  1 

3a3y  3p3a  3a  3p  3a3g 


32Q  32Q  .  3Ek  3Ek  32Ek 

-  «•  -  =  E[(  -  )(  -  )  +  Ek  -  ] 

3y3jj  3p3y  3y  3p  3y3p 
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